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Improved method and device for local spectral analysis of an ultrasonic signal 
DESCRIPTION 

Technical field 

The present invention relates to a method to process ultrasonic signals, 
5 especially in the field of non-invasive clinical tests based on the ultrasonic 
method. The invention also relates to an ultrasound device for the 
implementation of a processing method of the ultrasonic signals. 
State of the art 

Ultrasound scanning is a non-inVasive investigation method for 
1 0 examining internal organs of the human and animal body and also for non- 
invasive examination of other types of structures, including non-living and 
inorganic structures. As is known, it is based on the production of a series of 
ultrasonic signals (forming an acoustic beam) which are directed into the 
patient's body and which excite an ultrasonic response signal, caused by the 

15 back reflection and/or scattering of sound waves by the diffusers and/or 
reflectors present within the region under observation. The probe for 
transmitting and receiving the ultrasonic signal sends the response signal to a 
receiving and processing system, after which the signal is transformed into a 
two-dimensional ultrasound image. 

20 Commercial ultrasound scanners supply images relating to the 

amplitude of the ultrasonic signal back-propagated by the various tissues or 
interfaces struck by the acoustic beam. 

The various tissues or the various interfaces encountered are detected 
on the basis of their echogenicity, determined by their acoustic impedance, 

25 which in turn is the product of the density and the acoustic velocity. 

The purpose of conventional ultrasound scanners is to reconstruct and 
maintain the spatial order of the tissue under examination, by means of the 
amplitudes detected in the time domain, which in turn provides the spatial 
location in depth. The conditioning performed on the signal is aimed at 

30 maximizing the signal to noise ratio and keeping the sensitivity of the 
investigation constant when the depth varies. 

When the acoustic wave is propagated in a medium and encounters a 
non-homogeneous region, some of the energy is back-propagated towards 
the transducer and some is propagated into the tissue. 
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the envelope operation carried out on tne raaio .. m 
^o Uuencv s ig na, and therefore prov.es an es«mate of the s, g na, 
2S amP 'T h d e e 'enve,ope operation is current* used in a„ uitrasonio imaging 

stands for brightness) for the resuiting ,mages, * amp Wude 
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transducer consisting of a signa. having a carrier frequency approximately 

il a frequency modulation (asymmetrical spectrum) w,t phase steps 
frequency pea k s)and possibly w«h a morphology due to nonpar d*tortion 
5 Inmemediumshowsnon-linearenects^ductionofharmomcsp^a. 

FigU re 1 in a conventional A-mode representation. It will be noted that 
SLLn on the frequency and phase of the original signa, has been 

^TTaliar, . the radio-frequency is eliminated by an envelope 
extraction operation, there is a loss of: 

information obtained from non-iinear interactions between u^ ound 
and biological tissues (local acoustic impedance wh,ch depends on the 
density of the medium, which in turn depends on me .ntensrty of the 
15 incident opsonic wave, re-emission on the harmomcs and 

subharmonics by bubbles present in the contrast medium, etc ); 
. infcn.ationobtainedf.m.ine^ , 
signal spectrum due to seieotive absorption, especa.hr where th,s ,s 
pthJpendent. and the scalloping effect (indentation of the shape o 
20 the amplitude spectrum present in signals which have been 

LscaLd by discontinuities, including those of the pent type, 
distributed periodically in space). 

Iventiona, uKrasound scanning systems used for medica, d.agnos 
are therefore based on the observation of an A-mode trace or a B-mode 
25 image, which show only the amp.itudes of the echo signal relating to the 
various tissue structures encountered. ,„,„*,„*» 
Consequents there are situations in Which the conventional methods 
are found to be limited since, owing to low sensitivity or low specfiorty, 
Z ?ossib,e to obtain a reliable diagnosis of pathology. However . 
30 frequency signal contains not only the information relating to amphtude, but 
30 ToTe Uency and phase information (related to the time de,ays wrth 
which the echoes due to individual diffusers reach the transducer). 

The latest signa, processing methods tend to extract all the parameter^ 
characterizing the signal, in particular what is known as the spectra, content. 
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Spectra, analysis applied direct* to the .die frequency sfcna, -n^r** 
information from ai. the characteristics of the echo s,gna, as * ,s returned 

from the transducer. 

A spectra, analysis a,gorithm for "reading" the s,gna, ,n the frequency 
5 d omarnma k es«poss, b ,etoprovidecharacteristicva,uesford m erent,at,on^ 

Te various types of signa, and consequent* of the bio,ogica, structures wrth 

^rtin Figure 4 again shows the radio frequency signa, , 

Ficure 1 while the diagrams (B) and (C) show the spectra, content of the 

fori portions of the trace in Figure 4A, to show how the structure 

lences in I fissues resufl in different types of ultrasonic imposes and 

consequent* distinctive frequencyvanations in the respective spectrum. .V. 

spectrl, contents of me two portions of signa, shown n F,gures ^ nd 4C 

elnd from 0 to 12 MHz. ,t wi,l be noted that «<*~£j£ 
approximate,ymesamecentra,frequency(thefre<,uencyatwh 1 chthepea k ,s 

fou nd) in the vicinity of 4 MHz (the signa, re,ates to scannmg w* an 
masonic probe having a nomine, frequency of 5 MHz) but ave a afferent 
frequency variation, since the first graph has contnbuflons (above 5 MHz), 
which are present to a ,esser extent in the other portion of the signa . 

,n practice, the ufirasonic signais backscattered from the 
microstructures, due to the deferences in the propagafion P«ereon 
.he face of the receiving transducer in a destructive or consfructv e wa* 

of the echo signa,, that is, the echogenicHy. dependmg on the *pe 
interference, is a,so infiuenced by the refiectiv^ of the fissues w^oh « *n. 
determined by the mechanica, nature, densfty and elasficty o the 
dents i the mi— re; for examp,e, a high content - ^agen 
increases ,oca, mechanica, impedance and thereby echogemcrty. Therefore. 
h e lowing contribute to the appearance of the classic ufirason.o ,mag . 
pa , agglomeration of the microstructures, which determines the specie 
oLceriL, and their mechanica, nature which determine^ o gen,^ 
The mean amptade variation in space therefore seems not to P rov,de 
bharacterizafion of the fissues as , is influenced by venous 
morphological factors which are formed randomly. 
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if intonation oan be extracted on the speckie in terms of spatial 
periods (spatia, frequency) then ft is possib.e also to acu.ro .nformahon 
reiated to any spatia, periodicity of the agglomerations under 

The radio frequency uftrasonic signal is influenced by the mfonnation 
5 produced by the speckle of the images, in terms of particular amplftude and 
phase modulations of the ultrasonic signal. 

The information present in the time domain are preserved n *e 
frequency domain and the influence of the speckle is translated ,n o the 
scalloping effect that consists in a variation of amplitude accord.ng to *. 
,0 graphlLpresentation indicated in Figure 4 D . .n fact, for ^examp, ^ he 
Lse of N diffuse, or reflector* (the distinction, which ,s not unequ^oca. ,s 
made in relation to the wavelength of the ultrasonic sfcnaD struck by he 
masonic beam, the spectml density of power, calculated accord.ng to the 
Fourier transform, of the echo signal received \X(.f may be represented as: 

'mar =m\p- +2 1 %7' ~ h)h \ 

where: 

H4 2 : spectral density of power of the signal transmitted 
ri: reflection coefficient of the i-th reflector 
li: distance between the first and the i-th reflector 
c: velocity of the ultrasounds. 
,n the case of only two reflectors separated by a distance equal to 1 , the 
expression takes a simpler form: 

\X(a>f = \S(a>ffc W + 2r x r 2 cos(2^/c)] 
T m *L = E is the distance in time between the beaks of the signals 



15 



20 



Leo 



25 back-propagated by the two reflectors. 

ft can be seen from the schematic formulae that the spectrum ,s 

weighted byacosine with period proportional to the distance of the reflectors; 

this determines the minimum and maximum values of the sbectra, pow 

density This effect is called scalloping, and is represented schematically n 
30 £T* for two reflectors, posftioned at variable distances (indicated wfth 

SO S1 S2)fromeachotherinthedirectionofpropagationof«heu,t ra son,c 
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in order to distinguish two otherwise indistinguishable reflectors, 
at disance smaller than the duration of the ^^e 
to evaiuate the distance between minimums in the spectra ™ ^ 
case of severai reflectors the spectrum of the initial s,gnal .- -^ted by 
tapped osciiiations having different period, The fonnulae aiso « do 

-=^=^-=^^ 

the elements for differentiation of the biological structures. 

.oio-freguency signal is a powerful instrument for the character, of 
tissues by means of the relative ultrasonic signal. 

Tn Lent years me development of new signal process.ng algonthms 
has mad" P lie to obtain further information on tissues, by inveshga ng 
!T*n 1 the ultrasonic beam in terms of spectra, variations mduced ,n 
.nteractonw, it nas b^n found that pathological stmctures can 

— r»— — — — — — 
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However, the methods developed up to the present time ^ have not 

where the speed of data output is a basic requirement. 

E p A- 341 003 describes a new method for spectra, analys.s of an 
u ,trasonic signa./which overcomes the Hmits and drawbacks of prior art and 
30 rextraL of the information product the specie of the ^ 
content of this previous publication is fu.iy incorporated ,n the present 
description and forms an integral part thereof. 
ohi ^ c a nrt Summary of the Invention 

yBi§ss ^ D rTZTrrihpd in EP-A-1.341.003, the object of the 

Starting from the method described in a 
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present invention is to ^ an imp— and e« * • *• 
possibles offered by tbis type of spectra. ana.ys,s, to obta.n greater 
more precise information on the structure under exam.nabon. 

- in EP A 1 341 003 is based, the invention descnbed here>n allows 

— — »-«»»~^— i 
— — rr sr^r. - — - - 
rsr— > «- — • » - ™ M ,ri^ 
-rr— r:»::rx. 
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radio freauency signal returned from the structure examined. 

"be apparent from the text hereunder, the .oca, estimators are 
re prestn ed V the coefficients of a sufiable decree, for example fourth or 
Id L ee interpolating pofynomia., which approximates the coeffic,en, 
« Zt Tfo'r the various frequency hands into which the digfcze and 
saml radio frequency signal was previous* subdivided. Aitematveiy, the 
oca —re may Insist in combinations of coefficients of vanous 
teeToTte interring polynomia,. The combined of coefficien^ 
various de^es can take place before or after having weighted each 
coefficient with the shape factor of the distribution histogram. 

interpolating polynomial are utilized instead of a single coefficient. 

^sIms.™ into bands, that is, decompose into bands or sub-ba d^ 
can take place with various criteria. Moreover, also the shape parameter of 
20 Z dlln histogram of the loca, estimators can be chosen ,n vanous 

the T 9a l"m:rent levels of decomposition both to obtain dmerent 
25 spatial resoiutions and to use various subdivisions of the overall spectrum 

T '"T— ns of the coefficients of the interpoiating 
polyn omia, "weighted" with parameters related to the nature - ~ — ^ 
distribution by means of a targeted study of the histograms through specific 
30 parameters such as: 

most populated classes 
- standard deviation 

kurtosis 
. higher degree moments 
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experimental tests performed. Derform ing a combination or 

According to a different aspect, instead of performing 

tinn of the iocal estimators and the shape parameters of the 
weighting operation of the local es inV ention two or 

motion histograms of said estimators, -~^^n -n be used as 
mor e .oca, es«mators can be combine ioca, estimators 

a combined local estimator. For exampie, several — each 
can be P—d. in which ^£^,£2; polynomial, 
matrix, by coefficients of vanous degrees r ^ can be 

, Subsequently, coefficients in homologous pos^ons of these ^ ^ 
combined, for example the coefficients of all o some ^ ^ 
added or multiplied together, in concius.n ,n . 
is obtained, each of which is composed of a funct.cn 
contained In the previously obtained matnces. 
5 raetdescnjrttonj^^ lMr1v understood from the description 

. rl of the spectral content in two portions of sad trace, 

! ° and zr 5 r: - - — — 

3,5(53 eA -ows an u«rasonic excitation signal sent by the transducer 
into the organ to be examined; 
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Ftaure 10 shows the distribution histograms of the vano 
ooeftiST ^ the fourth degree interpoiating poiynontia, — - * - 

-CTl — Lograms of the ooe ffl cients of »e 

interpoiatingpoiynomiaiiothecaicification region g 
Figures 13A-13C show the venation of the ' «* 

Bgures14.15and16show,an,^ouatoBgures10^ 

. • x nf the various coefficients ot me 

distribution histograms of the gnd 
interpolating polynomial in the various reg.ons of hp* P*. 
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calcification; 
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calcification; distribution histograms of the 

bands or sub-bands by means of a discrete wav 

Figure 18B schematic^ shows the six decompose bands. 
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Figure 18C shows the variation of the four* order inflating 
poiynomia, obtained by applying the decomposition ,n F,gures 1 8A, 1 8B 
ultrasonic response signal of a carotid artery; and 

X es 19, 20 and 21 each show a B-mode uitrasound .mage of a 

nrocedure of the local estimators. 

L^om^?^^ uKrasound 

^Tti^aW^ch— "lock d-agram of a ufirasound 

nnaratus The number 1 indicates a transducer of the scanmng 

and receives the uftrasonic response signal. The analog output 3 of the 
ICead is connected to an acquisition card l The = produ- 
sampled and digitized radio frequency signal, wh.ch forms the ,nput for 
baTof integrated FIR fitters, indicated generic* with 7. The bank of fitters 
S I dn sist of any commercially available integrated device, for example 
7 can consul y duced by 

the 43168 circuit produced by Harns, 

^^T™ - ^ a controller for the £ 
processing unit The microcontroller 9 can be. for example, a TMS 
nSP Droduced by Texas (USA) or another equivalent. 

Th number 1 1 indicates a memory associated wtth the filter 7 and wtth 
the confer 9. The memory 11 is connected to an analog/dig^ converter 
13 and to the digital stages of the scanning head 2. 

The circuit described here is known and does not requ.re more detailed 
aescnSon Essentia,*, the apparatus operates as follows. The transduce 1 
a series of ultrasonic excttation signals having, for example, the form 
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• Fi „ 6A into the body containing the organ to be examined, for 

5 S ' 9na L radio freguency signa, is taKen from the receivin £ , cimurt of the 
scanner immediate, after the - 

and is then dig.* - -^^ti » - 
and with a resolut.on of 12 bits (or o 
,0 acg^on card , in P-ce ~~ 
until a frame or a porton of a frame i bank of the 

~- - - 

acquisition card until a compie processor 

25 of the signal spectrum. operations carried 

c- , mc 7A 7E summarize, in the form of graphs, the operauon 

rsrrJLr— rr — 

EP . A -t 341.003, where not described in detail here,n. ^„ 
Processing of the samp.ed and digitized signa! acqu.red by the card 5 
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sianal into sub-bands. Essentially, for each track ef the acquired, sampied 
7 i^ed Irne, the s,na, is decomposed into a series of frequency 
ba nds or sub-bands using a fime-frequenc, transform euch a a wav^e 
in oarticular by means of discrete wavelet packet 

bands or sub-bands, which wiU be indicated below as band 0. band 1. band 
7 Howeler, as shai, be apparent hereunder, this va!ue is not b.nd.ng; on he 

by filtering the digitized and sampled RF signal ,n F.gure , 6B ^However, the 
of structure or fis ^ ^ ^ . 

Has I —d w«h two sections in cascade. The coefficients of *. 
C which provide a wavelet transform are defined — ^ 
Lwn For exampie, « is possible to use the coefficents defined ,n US 
SZg 8 the content of which is incorporated in the present descnphon 

Zl 7A schematically shows a frame acquired by the scann.ng head 
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r -rh* frame referred to below as the "input frame", 

dis tribu«on in bands or sub-bands .<*«-. * ^ 

each of the bands shown schematic* >n F,gure 7C and .nd.cated 
15 band 1, ... band 7. f ^ same 

grea ter number o, decoding sub-bands from ^ o free J ^ 
grea ter amount of information on its spec*. ^ ^ 

bands can also be obtained w. H suW3an ds with a 

ampmudes) which are no. constant that ,s^ba ds o ^ 
variable amplKude and with filters of any form The *~ 
decomposition will follow experimental cntena as a functor 

^rirroTrrcr^sformia 

The matrix— 3*-^"^ coefficients of said matrix are 
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the dimens,on of t,me (that ® of d P , ^ 
0 _ 7 are taken along each *"* fe determined . This 

which approximates the variaton of the- coeftaents 
operation is shown schematic^ ,n F,gures 7D, 7E More p ^ 
70 shows the spectra, coefficients in the venous „ ^ 

, and to tracK n, by way of examp, tf a ^ aligned 

coefficients of the various hands *7 ,n pos* U and^ 
with one another respectiveiy in the first and ntheeoon 

7E , where Ph indicates ^^J^^ attnelpth „ 

:rr:r P p :^r — — • - 

Oneormorecoeff^o t J y ^ ^ ^ 

which allows, as will be explained D ^ 
inf orma«on on the spectrum of ^ ^^^^ of the 

20 perameter is ^oafive. the — n - ^ ^ ^ , ^ 

signal. This is because rt ■ related, hy coefficients in 

• p 0ly nom,a. which approximates the vanahon o the s^ec ^ 

given point of the input frame. he 
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structure to be investigated. Some non-limiting examp.es w,« be deson 



below. 
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The .east squares method of approximation is usefu. for the oaicuiation 

- cetermine the poiynomia, of degree , m <jr V « ^ 
experimental observations yo,y< , y„-i correspon 

X1 Xn-1,: m 

P(x) = a 0 + aiX+ + a m x 

approximate and the observed values, 

S-g(*b+«i*/ + +a^7-y/) 2 

the coefficients must be sited in sucn a way that S is minima ,,in order to 

-^=====^---- 

Si= f, ? and thesystemofe q uationsisobtained,whichoan 
be written in matrix form as: 

^t^^^^^T^™* for various 
It was observed for the various types of tesue 
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- —7 — t — Hn. that * *• 
distribution h.stogram of each coe pathological 

o, the histogram in various conditions and for vanous tissues, refer 
should be made to Figures 8A to 17B. . and 

ooefHcients of the poiynom,a. have « s,gnff _ 

20 ^r:rrr:r S orc:t — « 

the image are three Endows centered on the , , = ^ ^ 
ofmevessa,,tha,is,ontheb,cod,andon*e^o^ 

of the matrix of toed estimators were calcuiate d , th «e J9 
«ssue structures differing from one another, show, hat *ese es 
have a statistical distribution which vanes accord.ng to M> 

iWeS ££ 9A 9B and 9C repent the variations of the three fourth order 
Porvnomials bating ^^^1^ ■£ 
deifications. In other terms, followng the fourth order 
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the portions of image relating to a lipid plague, the vanat,en ,n F,g e 9BJc, 
the portions of image reiating to the free region of the vessei, ,n 
flows, and the variation in Figure 9C for the portions of rmage reiaflng to the 
reaions in which calcification is present. _ 

als C ns of a porflon of oarotid artery. The olasses of reference are 

For each of the Ave coeffloients taken into consrderaton (mdrcated hem 
wrth AO A1 A2 A3, A4), three histograms were oonstruoted correspond^ to 
II isfl L — ositioned respective, in a region of the lipid plague 

order coeffioen, (A2) jn papula ^ ^ ^ ^ 

distribution compared with that ot me be»u« 

a d calciflcaflon. In other terms, hy anaiyzing the 
h t nrams representing the statistical distribution of the second order 
" Zoll distinguish whether the region examined contains 
Vessei (blood, Figure 11), lipid piague (flgure 10) or calcifloabon (Frgure 

12> ' This means that iipid plague and caloiflcaflon introduce d«ferer,t 
spectra, modifications in the back-propagated signal, which can be detected 
J^lave evaiuaflon of the distribuflon histogram of the second orte 
events of the Interpol polynomia. which approximates the vanatron 

^ the r 9 r — ons can be made regarding the variaflor . of the 
■ «nn to the third order coefficient, which has a different 
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results a h the varia tions of the fifth order 

coefficients. Figures 13A, 13B and sn 

distribution of the coefficients (A0-A5) ot 

, fnr the same windows indicated in Figure ba, 

interpolating polynomial for the sam on 
centered respectively: on the lipid plaque, on the blood 

affected by calcifications. djfferent 

nl*nd The B-mode ultrasound image used to obtain the 

h »»<th V tissue and neoplastic tissue, also for the shape of the 
between healthy fissue P popu ,ation frequency 

distribution histograms. Figures 17A and 17B show the p 
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rFiaure 17A) and for the coefficients a3 and a^ ^ny 

healthv tissue, and moreover in the case ui K « eQOn 

L t» °f »■ **" ™ 
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the matrix of the local — - - ~ 
vanousordersofagrven p y ^ jmage examlned 

tocal est,mators. a senes o *"* indjcating wfth „ 

equa. to the number of coefficients a, of the matnx of .oca. est,mato 
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value of which has been modified and prooessed taking into con— 
r hape palter (in the specific example the "standard dev,a on * of 
ZZL histogram. This matrXwil, - composed of m x n values b 9 
bij = f(aii,og) 

* Qr ^„ For example, each coefficient ay of 

eg obtaining: 

■ rrrr:rr.rs= 

plurality of matnces of local est,ma distribution 
weighting them with the shape parameters of the relet, 

hiSt Th7shape factor or shape parameter of the histogram which is used 
to use the kurtosis, the symmetry .ndex, h.gher degree mo 

here u2 c2 matrix of weighted loca, eafimators. 
^cient of said matrix has been "weighted by a she e ^ 
30 distribufion histogram of the ^^^^Z^ 
. — "-"TeLon or another suable 
parameter. 
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„ must be understood that the operation described may be performed 
separly for al, or for some of the coefficients of the various orders o *e 
ZTJng pofynomia, which are usefu, for speotra. charactenza on o he 

of the various orders of interpolating polynom,al ,n drfferent ways. For 
L^e the third order coefficients oouid be weighted by <« by ffi 
Zeter of standard devia«on, whi.e the second order £J 
efther be weighted, or another operation of oomb,nat,on wrth a dftferen 
rietlr, such as Kurtosis, or an operation differing from dtvid.ng the ,oca, 

account the shape of the distribution histogram of sa,d coefficients 

SPe *gie schema indicates the procedure which can be 
perJmed when, instead of a singie coefficient of .a ,n.« 
20 porynomiai. aii the coefficients of said polynomiai are used to obta n a ^more 

rsdmatorsistherefor.amree^imensionalmatrbcwnhd.mens.os. n 

digWz ed and sampied uftrasound frame and K is the number .f c— of 
25 the interpolating polynomial (polynomial of degree k 1). 
the matrix are indicated with aii (1) , ...ay . 

tnemaxrx com bined with one another to 

The coefficients a^\ a„ , ...a,, are 

=" 

^ f ge neric function. For exampie. ft may be a summation or a product of 
ailthekcoefficientsofpositioniiofthethree-dimensionalmatnx 

On this matrix c, the shape coefficients o> for exampie the standard 
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The proceau France that coefficients (weighted 

_ _. mts 7n 7 u althouqh with the ditterence ui«i ^ 

the interpolating polynomial, instead of on* one of these, 

10 final matrix. coefficients of the 

intermediate solutions, in which only some of the coeffic, 

' weighted local estimators are then ^collated 

bij^ = f(ag W , aft 
Finally, the weighted localestimators are combined to obtain 
c i = f(b, <1> b,«) 

■--rr.s 

,onger processed to obtain the shape factor o* but its coeff.ce 
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estimators, not weighted) are suited to stafisfica. an^is to direcfiy 
obtainthepopulafionmatrixjndicatedherewithf,' es 
Having obtained the matrbc of the weighted iooal estimators, the values 
of the coefficients of this matrix are Casstfed according o 

• ««+oiu/ in be reated bi-univocaliy to a 
5 which have been *und * proven that a 

specific tissue structure. In other terms, P Droste te) 

certain type of tissue (for example the neoplasfic fissue of the p»*W 

points of the matrix of the weighted loca, estimators, fine value MM 
within the characteristic class of the neoplastic tissue, ,s overla.d on tine black 
2 white B-mode image. This color indication provides the , *. nd 
operl wfih the required information, showing the presence and position of 
"7; — — — - -perfonmedforvadous^esof 

^^rrrvrr.^ 

iloe the presence of two different types of tissue or structure (for 
same image the presence u pnouah to 

example lipid plaque and calcification in the carotid artery), rt .s enough 

rrr— " 
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•™ «f h^althv tissue, adenocarcinoma and 
example, the homogeneous regrons £ h hV ^ ^ 
adenofibroma are identffied. The rad.o frequency 
in response from soanning of the experimental samp.e, as ^descnbe i a 

- i rrs=r— — ~:r, 
„.„ » «— - *~ ^^jz—. . ~ 

procedure will be used for the otner iy P 
adenofibroma.inthecaseoftheprostate). 

More detai,ed processing envisages wodun on atleas, P 
ma tdces of local abators both in the ^ound 

Acguis.cn of an uitrasound frame ,s ^ X JeT l to obtain 
experimental sample, which is then processed as descnbed 

looal abators, for exampie those ^ - ^ 
coefficients of the interpolating polynomial. A second g P 

^^^r^niSTli. mat. of 

- - o.er efficients of »e 
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■ , w =, nrrctate The coefficients of the matrix of 
interpolating polynomial for a prostate, me c 

weighted loca. estimators obtained by order 0 coefficents of the same 
Zl. are shown on the ordinate, Figure 23B shows, on the abscissas 
the weighted loca, estimators obtained by second order coefficents and, on 
the ordinates, those obtained by order 0 coefficients. 

Therefore, a pair of values (the coordinates of the pent), compose by 
*. value of one and of the other of the weighted .ocal esttnators used, ,s 
associatedwitheachpointofeachofthetwographs. 

The graphs obtained in this way can be overiaid on the >npu, unsound 
franJ ThereL it is possible to associate with each point of the rnput 
1 pa. I. oiasses of values of the weighted local estimators 

examined and corresponding to this point. 

By means of histological analysis, the regions of the tissue wh,ch « 
h0 mog neous in nature (healthy, neop.as«c, etc.) are ^ed on the 
Zs und frame. In Figures 23A and 23B these regions are «^ 
three different shades of grey. As a pair o, values of the M~gM- MocaJ 
^meters correspond, in both graphs, to each point of the ultrasound frame, 

1 each homogeneous region of the tissue examined. sa,d class be.ng 
20 univocallyassociatedwiththatspecifictypeoftesue. ^ 

,n order to identtfy each region bi-univoca„y, it is important to select^ 
set composed of pairs or tripiets of classes of values of the we,ght« 
elators, that are not overlaid on any other associated wrth a drfferent 
; e ;i The operation to extract these sets can be performed through, for 

25 — T S^tSX- of *e sets thus iden.ed is 
represented with a color scale, dtferent for each set, on the conventional B- 
leCe and * direct* proposal to the intensity of the population 

30 ""ad mentaliy, as can aiso be noted in Figures 17A, 17B, the 

which is desired to show (in the exampie the neoplastic tissue can non,a y 
I highiighted by combining simp* with a sum the images obtained from 
the coefficients of the interpolating polynomial which can be h,ghl,ghted by 
in this case, w*h an AND opera«on the sets of values obtained 
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iJZZL. detected * anafysis and Casses of va = 

o 1 weighted local e— > aiso va,id in the case 
5 ioca. estimator for example obtained by combing deferent degree 
coefficients of the interpolating polynomials with one another. 

T e experimental resoits discussed above and indicated ,n £M- 

^^^^^^ 

le^n Tl prostate gland typicaify has structures with dimensions 
SZ, - carotid adery. The deferent dimension of the structures , that 
"fleet the opsonic exeflation signal cause differeht distnbubon of the 
Z^L Ze return signal. Subdivision into sfcteen bands of the same 

RF siqnal returned from both the organs examined. 

Nonetheless, in genera, i, wil! be understood from tte 
that fine ohoiceofdecomposition into bands(number and 
may differ according to the type of tissue or structure to be ana^The 
smaller the structures to be studied in the analyzed organ are. the greater the 
Jsoln of decomposition into high frequencies will be, and v,ce versa. Tn* 
r^Ld by the experimental results discussed below wflh reference to 
B^Ts to 21, obtained again from processing an ufirasound ,mage of a 

P °XZTT^ shew the B-mode image of the podion of caret, 
e^lyzed. The digged and sampled BF - < ^ = 
discrete wavelet into stx bands of variable amplitude. F,gure 18A shows t 
te representing the filter process used and Figure 18B ; show^ the « 
• + - h.nrt. It can be observed that the amplitude of the bands 

IClg PoVnomia, obtained applying the aforesaid decompose to the 
B-mode image in Figures 19, 20 and 21. 
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Figure 1 9 shows, in the top part, the B-mode image of the caro d arte* 

Responding to Figure 18C) and subsequently the distribute h,stogram 
A0-A4 of the various degrees of interposing poiynomra 
ZZ at the ieve, of a window (shown overiaid on the B^e = , 
centered on the lipid plaque. Figure 20, with the same layout of the ■mages, 
rhlTheTntenooL; polynomla, and the distribution histograms « 
awindowoentered on me free region of me vessel, that ,s, on thebood^so 

ease the window is shown overlaid on the B-mode image of the carofid 

rriSLL- on a portion of the B-mode image in which calcificafion 
phenomena are recognized. 

* a on and 21 t can be seen tnai wiui 
By comparing Figures 19, 20 and zi 

frequencies, the variation of the interpreting polynom,al does not change 
ZLw in the three regions observed (lipid plaque, biood and 

roeTnot provide resolution appropriate for examination of this type of 
2 1 due to the fact that .he dimensions of the structures to be recogn^d 
20 Ismail and, therefore, the spectra, information they supply is located at the 
" ^equencies of the response spectrum, that is, where the ^eoompo^on 
used requires bands with a greater amplitude. Contra* to *h,s the res Its 
obled on the same organ with decomposition into sixteen bands o he 
sllpmude (Figures 9A to 16) show that the resolution obta,ned w,th th« 
decomposition is sufficient to provide significant results. 

7L the distribution histograms in Figures 19, 20 and 21 t can be 
noted that in fact (although the variation of the interposing polynom,a. . 
ienfially the same in the three regions of lipid plaque, blood and 
the shape of the disfribufion histograms of the various orders of 
30 ZSL changes. Nonetheless, the deference between d,stnbuhon 
^remsofhomoiogcuscoefficientsinthethreeinvesfigatedreg.ons.s o 

a marked and significant as the difference obtained using sighal processmg 
Zfon decoding into sfcleen bands of idenfica, amplitude, n conc,us,o . 
caparison of me results obtained with the two types of decod.ng into sub- 
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Figure 22 snows Transform 

H *, le characteristics of the examined structure, the investigator, method 
Z 1 an" the various decompose P oss«ies a,on g the tree ,n 

« the drawing shows on* one practice, embodiment of the 

any reference numhers in the attached a* - «*p 

facilitating reading with reference to the P reo f 109 
drawings, and does not iimft the scope of protect™ of the cia,ms. 
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